Phase behavior and material properties of hoUow nanoparticles 

U. S. Schwarz and S. A. Safran 
Department of Materials and Interfaces, Weizmann Institute of Science, Rehovot 76100, Israel 



O' 

o 

(N; 

Oh- 



O 



O 

o 



> 

o 
o\ 
o 
o 
o 

-(— • 



o 



% 



Effective pair potentials for hollow nanoparticles like the 
ones made from carbon (fuUerenes) or metal dichalcogenides 
(inorganic fuUerenes) consist of a hard core repulsion and a 
deep, but short-ranged, van der Waals attraction. We investi- 
gate them for single- and multi- walled nanoparticles and show 
that in both cases, in the limit of large radii the interaction 
range scales inversely with the radius, R, while the well depth 
scales linearly with R. We predict the values of the radius R 
and the wall thickness h at which the gas-liquid coexistence 
disappears from the phase diagram. We also discuss unusual 
material properties of the solid, which include a large heat of 
sublimation and a small surface energy. 



I. INTRODUCTION 

The phase behavior and material properties of atomic 
as well as colloidal systems can often be understood using 
the concept of an effective pair potential |p,p|. However, 
in the presence of attractive interactions there is a fun- 
damental difference between the atomic and the colloidal 
case: while atomic systems usually have a range of at- 
traction which is comparable with or larger than the hard 
core diameter, colloidal systems have a range of attrac- 
tion that can be much smaller. One example is the attrac- 
tive depletion interaction in systems with sterically stabi- 
lized colloids and non-adsorbing polymers. The range of 
interaction can be tuned by the polymer's radius of gyra- 
tion and usually is one order of magnitude smaller than 
the hard core diameter [p|. Another example is the na- 
ture of attractive surfactant interactions between inverse 
microemulsion droplets. Here the range of interaction is 
fixed by some microscopic interpenetration length. By 
tuning the droplet size, that is the amphiphile/oil ratio, 
the hard core diameter can be made more than one order 
of magnitude larger than the interaction range ^ . 

For Lennard- Jones systems, the potential depth de- 
termines the temperature and the hard core diameter 
the density scale, respectively; after rescaling to reduced 
units for temperature and density, no other degrees of 
freedom are left and phase transitions fall on universal 
curves (law of corresponding states). The resulting phase 
diagram features a critical point for gas-liquid coexis- 
tence and a fluid-solid coexistence which is first order at 
all temperatures. For more complicated interaction po- 
tentials, the topology of the phase diagram can change. 
It is long known that the gas-liquid coexistence disap- 
pears for small interaction ranges, as is the case for the 
short-ranged depletion interaction in systems with steri- 
cally stabilized colloids and non-adsorbing polymers [0. 
Recent computational M and analytical MM work sug- 



gests that the liquid phase disappears and there is only 
vapor-solid coexistence when the range of attraction is 
smaller than about one-third of the hard core diameter. 
Moreover, an isostructural solid-solid transition appears 
when the range of attraction decreases by another order 
of magnitude. 

Recently, the phase behavior of Cgo (buckyballs) has 
attracted much attention since it constitutes a border- 
line case for the disappearance of the liquid phase. The 
van der Waals (vdW) interaction between buckyballs is 
usually modeled by the Girifalco potential which is ob- 
tained by integrating a Lennard- Jones potential for car- 
bon atoms over two spherical shells ||] . It is expected to 
be valid above 249 K, since then the buckyballs rotate 
freely; for lower temperatures, the interaction becomes 
anisotropic Q. Since the length scale for attraction is 
set by the position of the minimum of the Lennard- Jones 
potential (3.9 A), it is smaller than the buckyballs' diam- 
eter (7.1 A). Early simulations using this potential gave 
conflicting results in regard to the existence of a liquid 
phase: while molecular dynamics simulations predicted 
that a liquid phase exists pOJ, a Gibbs ensemble Monte 
Carlo simulation predicted that Ceo might be the first 
non-colloidal substance found which docs not have a liq- 
uid phase at all |1^. Recent Monte Carlo simulations 
however confirm the existence of a small liquid phase re- 
gion in the phase diagram P4[l3[ . 

However, carbon buckyballs are only one example of 
hollow nanoparticles. Until now, more than 30 other ma- 
terials which are similarly layered have already been pre- 
pared as hollow nanoparticles (either spherical or cylin- 
drical), including e.g. the metal dichalcogenides MeX2 
(Me = W, Se, X = S, Se), BN, GaAs and CdSe @. 
In fact, the formation of closed structures is generic for 
anisotropic layered materials of finite size due to the 
line tension resulting from dangling bonds. Effective 
pair potentials for such hollow nanoparticles are isotropic 
for spherical shapes, but still depend on their radius 
R and thickness h. In particularly, for carbon onions 
and hollow metal dichalcogenides nanoparticles (inor- 
ganic fuUerenes) the thickness h can vary because the 
particles are multi-walled. Carbon onions with hundreds 
of shells have been observed 111] . Onions are formed by 
metal dichalcogenides with up to 20 shells m&. In both 
cases, outer radii R can reach 100 nm, that is several or- 
ders of magnitude more than buckyballs with R — 3.55 A. 
Unfortunately, control of size and shape is still rather 
difficult, and only some fuUerenes can be produced in 
a monodisperse manner and in macroscopic quantities. 
The best investigated case are of course buckyballs, which 
readily crystallize into fullerite. However, the synthesis 



of hollow nanoparticles is a rapidly developing field, and 
it is quite conceivable that more experimental systems 
of this kind will be available soon, since hollow nanopar- 
ticles provide exciting prospects for future applications 
in nanoclectronics and -optics, for storage and delivery 
systems, for atomic force microscopy and tribological ap- 
plications p4[ . 

Although this paper is concerned with hollow nanopar- 
ticles made from anisotropic layered material, there are 
other types of hollow nanoparticles which should show 
similar physical properties. Colloidal core-shell particles 
with a porous shell can be made hollow by removing the 
core through calcination or exposure to suitable solvent. 
Hollow polyelectrolyte shells have been prepared by col- 
loidal templating and layer-by-layer deposition ^^ . The 
same technique can be used also to synthesize hollow sil- 
ica particles ||l^. Since wet chemistry and electrostatic 
self-assembly allow for good control of size and shape, 
the wall thicknesses h could be tuned from tens to hun- 
dreds of nanometer by varying the number of deposi- 
tion cycles; the radius R is determined by the size of the 
templating polymer particles (usually in the sub-micron 
range) . Hollow nanoparticles also occur in biological sys- 
tems. Clathrin coats are hollow protein cages which in- 
duce transport of specific membrane receptors into the 
cytoskeleton by coating budding vesicles, but which also 
self-assemble in a test tube. Recently, the highly regular 
structure of the smallest clathrin coat (radius 35 nm) has 
been resolved in great detail by cryo-electron microscopy 

In this work, we address the question of how the phase 
behaviour of hollow nanoparticles depend on radius R 
and wall thickness h. For this purpose, we use the same 
framework as the work based on the Girifalco potential, 
that consists of integrating a Lennard-Jones potential for 
single atoms over the appropriate geometries and inves- 
tigating the phase behavior resulting from the effective 
potential. We also predict how the material properties 
of the solid vary with R and h and discuss how elastic 
deformations will modify our results. 

Our main result is that for both single- and multi- 
walled nanoparticles, in the limit of large radii, the inter- 
action range scales inversely with the radius, i?, while the 
well depth scales linearly with R. For the phase behavior 
this means that the gas-liquid coexistence will disappear 
with increasing R. Our full analysis predicts that this will 
happen for single- walled nanoparticles around Re = 12 A 
and for completely filled nanoparticles around Re = 35 A. 
It follows that buckyballs, which have R = 3.55 A, are 
far from losing their coexistence region. We find that 
the effect of the wall thickness h is rather small: Re 
initially increases with increasing h, but then levels off 
on the atomic scale set by the Lennard-Jones potential 
between single atoms. Moreover the temperature range 
is essentially set by R and is hardly affected by h. All 
our results are derived twice, once from the analytically 
accessible Derjaguin approximation and once from a nu- 
merical treatment of the full potential. We find that the 



Derjaguin approximation predicts the right trends and 
provides physical insight into the underlying mechanisms, 
but overestimates both potential depth and interaction 
range. 

We also show that crystals of hollow nanoparticles will 
have unusual material properties. Their heat of sublima- 
tion (which is already unusually high for fuUerite) scales 
linearly with R; at the same time, their surface energy 
(which for fuUerite is already smaller than for graphite) 
scales inversely with R. Thus with increasing R it gets 
more and more difficult to melt the crystal, although 
it becomes more and more unstable in regard to cleav- 
age. The thickness h does not have any significant ef- 
fect here since the particles are nearly close packed and 
the vdW interaction saturates quickly over atomic dis- 
tances under adhesion conditions. However, we show 
that an important effect of the thickness h is to sup- 
press elastic deformations for increasing h. If one wants 
to exploit the unusual propcrites of hollow nanoparticles, 
multi-walled variants are favorable as they suppress elas- 
tic effects which will prevent gas-liquid coexistence and 
crystallization. 



The paper is organized as follows: in Sec. HA we dis- 



cuss the Lennard-Jones potential for the interaction be- 
tween single atoms. In Sec. II B we integrate it over the 



geometries of single- and multi-walled nanoparticles to 
find effective interaction potentials, which differ consid- 
erably from the initial Lennard-Jones potential due to 
the additional length scales introduced. Analytical pre- 
dictions for the potential depth and interaction range can 
be found using Derjaguin approximations, which describe 
well the relevant part of the full potential in the limit of 
large R and which are derived in Sec. II C . In Sec. Ill the 



Derjaguin approximations are used to predict the scaling 
of the material properties of the crystal with radius R 
and thickness h. In Sec. IV we investigate both the full 
potentials and their Derjaguin approximations to predict 
at which values of R and h the gas-liquid coexistence 
will disappear from the phase diagram. In Sec. |y| we dis- 
cuss briefly the effect of elastic deformations, and finally 
present our conclusions in Sec. VI. 



II. EFFECTIVE INTERACTION POTENTIALS 

A. Lennard-Jones potential for single atoms 

The calculation of vdW interactions between macro- 
scopic regions of condensed matter is a well-investigated 
subject in colloidal science ||2O,0,pl|,p2l . It is well known 
that the geometrical aspect of this problem is well treated 
by a pairwise summation of the microscopic interaction; 
the strength of the interaction {Hamaker constant) can 
be derived from Lifshitz theory or from comparision with 
experiment. The 1/r^-vdW interaction and its crossover 
to Born repulsion at small separation is usually well de- 



scribed by a Lennard- Jones potential: 
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Girifalco obtained the effective interaction potential be- 
tween two buckyballs by integrating Eq. (m) over two 
spherical shells ||^. By fitting resulting predictions for 
energy of sublimation and lattice constant for fuUerite 
to the experimental results, he found the values A = 
32 X 10"*^° erg cm^ and B = 55.77 x 10-1°^ erg cm^^ 
for the Lennard- Jones interaction of carbon atoms. Very 
similar values have been extracted from an analogous 
procedure for graphite sheets ||2^. For the following it is 
useful to characterize the Lennard-Jones potential by its 
hard core diameter a = {B/Ay/^ and its potential depth 
e = A'^/4B. For carbon, a = 3.47 A and e = 4.59 x 10"^^ 
erg = —0.11 kT (where k is the Bolzmann constant and 
T room temperature). 

In order to use the atom-atom potential from Eq. (U^) 
within the framework of a continuum approach, one has 
to know the density of atoms. For single- and multi- 
walled nanoparticles, area density g and volume density 
p will be used, respectively. The in-plane elasticity of a 
carbon sheet is so strong that its in-plane structure is 
hardly changed when the sheet is wrapped onto a non- 
planar geometry. For example, the area density of carbon 
atoms in a buckyball, g — GQ/AttR^ — 0.38 A~^ with ra- 
dius R = 3.55 A, is very close to the one in graphite. The 
area density for both WS2 and MbS2 is g — 0.35 A~^; it 
is close to the value for carbon, since the effect of larger 
bond length is offset by the existence of triple layers. For 
future purpose we define a dimensionless area density 
of atoms, r = ga'^ 
of atoms, X = P'^" 
X — Ta/l, where / is the distance between layers. For 
carbon and inorganic fullerenes, I = 3.4 A and I — 6.2 A, 
respectively. For carbon one then finds r = 4.56 and 
X = 4.60. For the metal dichalcogenides, r is roughly the 
same, but x is smaller by a factor 2 when compared with 
carbon. In any case, the dimensionless densities squared 
roughly amount to one order of magnitude. 

Although the molecular structure of the metal 
dichalcogenides is quite different from that of carbon, 
the values for a and e are expected to be similar. From 
the Lennard-Jones potential given above, the effective 
Hamaker constant for carbon follows as ir^p^A = 3.9 x 
lO^^'^ erg w 100 kT. Indeed this is the right order of 
magnitude for the Hamaker constant of vdW solids in 
vacuum [ pl[ , thus in the following we assume that the 
values of a and e given for carbon will give the right or- 
der of magnitude results for the metal dichalcogenides, 
too. 



B. Full interaction potentials 

When integrating Eq. (|l|) over the volumes of two hol- 
low nanoparticles, we introduce two additional length 



and a dimensionless volume density 
^^ The two densities are related by 



scales: the radius R and the thickness h. Fig. |l| depicts 
the geometry of two multi-walled hollow nanoparticles. 
For the following it is useful to define the two dimen- 
sionless quantities rj = 2R/ a and v = h/ a, that are the 
particles' diameter and thickness, respectively, in units 
of the Lennard-Jones hard core diameter. Since h < R, 
we have 2iy < rj. The Girifalco potential for two single- 
walled nanoparticles of radius R separated by a distance 
r follows from integrating Eq. (nh over two spherical shells 
g. Using the definitions given above, it can be written 
as 
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where s — r/2R is the distance in units of the parti- 
cles' diameter. Here i' — h/a does not appear since 
the shell is assumed to be infinitely thin. For buck- 
yballs, 77 = 2R/a = 2.05. Then the potential min- 
imum is at To — 10.06 A with a potential depth of 
Vo = -4.44 X 10-13 gj.g = -10.73 kT. In Fig. |a we 
plot the Lennard-Jones potential from Eq. (|l|) for carbon 
atoms and in Fig. gb the Girifalco potential from Eq. (g) 
for buckyballs. The integration over the given geome- 
try has strongly changed the character of the interaction 
potential: it is now two orders of magnitude stronger 
and rather short ranged. The concept of a small interac- 
tion range will be quantified below. The small width of 
the potential well becomes evident when compared with 
a Lennard-Jones potential with same effective diameter 
and potential depth (dashed line in Fig. ^p). The geo- 
metrical effect on the effective interaction potential is a 
subject well-known from colloid science. 

We now calculate the vdW interaction between multi- 
walled hollow nanoparticles, that is between two thick 
shells. We consider two balls i?i (i = 1, 2), each of which 
consists of a thick shell Si and a core Ci, which is a 
ball with smaller radius. Then the interaction between 
the two shells can be expressed in terms of interactions 
between several balls: 

*^Bi,B2 ^Ci,C2 
- iVci,B2 - Vbi.Ca) - (^Bi:C2 - ^Ci,C2) (3) 

'-61,-62 ' Ci ,0*2 Ci,-B2 '-Bi,02 ■ 

For two identical thick shells, the last two terms are iden- 
tical. For two identical shells of radius R and thickness 
h, we can write Eq. (H) as 

Vpir) = VRRir) + VB-h,R-hir) - 2VR,R^h{r) . (4) 

In order to proceed, we must integrate the Lennard-Jones 
potential from Eq. (H^) over two balls of unequal radii. 
The two integrations can be done analytically, but lead to 



rath er le ngthy forr nula e which are given in Appendix H. 
Eq. (|A^ ) and Eq. ([A3| ) together then yield 



which can be used in Eq. (||) to obtain the full potential 
Vpir) in analytical form. It diverges dX r — 2R, has a 
minimum at intermediate distances and decays at large 
distances as 



VAr) = -7:\\^-^^H4r.' - Qu^ + 3^/^)^ 
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which is the vdW interaction between two bodies of vol- 
ume 47r[i?^ — [R— ft-)^]/3 each. In the two-dimensional 
limit /i ^ 0, the full potential Vp between two hollow 
nanoparticles becomes the Girifalco potential Vq from 
Eq. (ph between two spherical shells (where the dimen- 
sionless volume density x relates to the dimensionless 
area density t by x = ra/h). In general, the full poten- 
tial Vp is very complicated; in order to gain some physical 
insight, it is useful to consider its Derjaguin approxima- 
tion. 



C. Derjaguin approximations 

The Derjaguin approximation relates the interaction 
between curved surfaces to the one between flat surfaces 
if the distance r — 2R between the curved surfaces is 
smaller than the radii of curvature R, i.e. if s <C 2. The 
interaction energy per unit area between two planar films 
a distance z apart and each of thickness h can be easily 
calculated to be 
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The Derjaguin approximation integrates over one of the 
two surfaces and evaluates W{z) for z being the distance 
to the nearest point on the other surfaces; for two sphere 
of equal radii R it reads 



VD{r)=TTRl dzW{z) 



(8) 



Using Eq. (|) in Eq. (|) we find 
1 / 1 
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We now discuss two limits of this result. For j/ <C 1, 
that is /i — > 0, we obtain the Derjaguin approximation 
for the Girifalco potential from Eq. (0): 
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where we have used x — Ta/h in order to get the two- 
dimensional limit. The same result follows formally by 
expanding each of the two terms of the full potential from 
Eq. (0) separately around s = 1. The potential minimum 
in the Derjaguin approximation is found by minimizing 
in Eq. (fol) for s: 
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The Derjaguin approximation is valid for s ^ 2; thus 
it will describe the potential well correctly if sq ^ 2, 
that is if ?7 3> 1. From Eq. ( pT| ) we see that in this 
limit of large R, the equilibrium distance of closest ap- 
proach, 2R{so - 1) = (2/5)i/V = 2.98 A, is indepen- 
dent of R. Moreover it follows that in this limit, the 
potential depth scales linearly in R. In Fig. ^ we plot 
the Girifalco potential and its Derjaguin approximation 
rescaled by Vo and for 77 = 2.05,4 and 10. While in the 
first case (buckyballs) the difference between full poten- 
tial and approximation is still considerable, for the last 
case it is already very good. Comparing values obtained 
numerically from the full potential with the values from 
its Derjaguin approximation from Eq. (pT|), we find that 
the approximation for sq is very good even for buckyballs 
with R = 3.55 A {rj = 2.05). The approximation for Vq 
deviates by 34 % for buckyballs, but only by 10 % for 
R = 13.74 A {tj ^ 7.92) and by 1 % for R = UTAA A 
(77 = 84.98). 

We now consider Eq. (||) in the limit 1/ 3> 1, that is /i ^ 
a. This is the Derjaguin approximation for a nanoparticle 
which is completely filled: 
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Equilibrium distance sq and potential depth Vo again fol- 
low from minimizing for s: 



If we compare this result with Eq. (|11|) for single-walled 
nanoparticles, we see that apart from the numerical pref- 
actors and the change from area to volume density, for 
single- and multi-walled nanoparticles the Derjaguin ap- 
proximation basically yields the same results for equilib- 
rium distance so and potential depth Vq ■ In particularly, 
in both cases the equilibrium distance is independent of 
R and on an atomic scale, and the potential depth scales 
linearly in R. The reason for this agreement is that in the 



situation of close approach described by the Derjaguin 
approximation, the size of the gap between the particles 
is approximately a. Since the vdW interaction between 
two films saturates on the same length scale, it doesn't re- 
ally matter how many walls there are. Moreover there is 
a cancellation of two effects: in the case of multi-walled 
particles, more matter is present and the vdW energy 
increases. However, in our description of thick walls, 
matter is now distributed uniformely in space and the 
effect of the first layer is smeared out, so the vdW energy 
decreases. In reality both descriptions are idealizations; 
however, since they essentially yield the same result for sq 
and Vo, the overall description of the potential well should 
be sufficient. Note that in both cases, the temperature 
scale Vo/k (where k is the Boltzmann constant) is ap- 
proximately SrjTji, where Tr denotes room temperature. 
Thus if a gas-liquid coexistence exists for the systems 
under consideration, it will do so at several thousand 
Kelvin. In fact the detailed treatments for buckyballs 
predict the gas-liquid coexistence to exist close to 2000 
K. 



III. MATERIAL PROPERTIES 

A solid forms at high densities or low temperatures 
since then close-packed particles can benefit from the 
attractive interaction without losing too much entropy. 
Many of the physical properties of solids are essentially 
determined by the properties of the potential well. We 
showed above that for large radius, R, the potential well 
is well described by the Derjaguin approximation. We 
now proceed to predict the material properties of such 
solids. By using the Derjaguin approximations rather 
than the full potential, we are able to derive analytical 
formulae which show the scaling with the different pa- 
rameters involved; although they are strictly valid only 
for 77 — 2R/<7 ^ 1, our results predict the right trends 
and orders of magnitude even for small radius. For more 
precise results for small R (relevant to buckyballs), it is 
easy to treat the full problem numerically ^. Because 
the Derjaguin approximation resulted in very similar re- 
sults for the potential well of single- and multi-walled 
nanoparticles, it is sufficient for the following to use the 
potential for the former case, i.e. Eq. (|ll|). 

Fullerite has several unusual material properties. It 
has a very high heat of sublimation (40.1 kcal/mol, 5 
times larger than typical for vdW solids) H] and it is 
the softest carbon structure; its volume compressibil- 
ity 7 X 10"^^ cm^/dyn is 3 and 40 times the values for 
graphite and diamond, respectively pq ]. We now investi- 
gate how these quantities scale for different R and h. The 
calculations proceeds as for Lennard- Jones solids p6[ . 
We assume the crystal to be face centered cubic (fee). 



The energy per particle is 

i 

-1-6 F(V2 r) -I- 24 y(%/3 r) -f 



(14) 



where the sum is over all fee Bravais lattice sites (except 
the origin which is occupied by the particle itself) and we 
grouped the lattice sites into shells {i} of equal distance 
{ri} to the origin, where shell i contains rii different sites. 
We explicitly give the nearest, next-nearest and next- 
next-nearest neighbor shells. Here r is measured in units 
of r„„ , the nearest neighbor distance in the crystal, which 
follows from du{r nn) / dr = 0. In the following we will use 
the nearest neighbor approximation, u w QV{r). Then 
r„„ — 2Rsq and u = GVq with sg and Vq from Eq. (|l] 
The lattice constant is then given by a = \/2rnn- 

The heat of sublimation per mol is simply AH 
NAu{rnn) = QNaVq with Avogadro's number Na 
6.02 X 10^3 jnol-i. Therefore 



(15) 



AH = ^X^Tr^NAT^T^ 



Since in the nearest neighbor approximation the heat of 
sublimation depends only on the direct interaction be- 
tween the particles and not on any geometrical aspect 
of the crystal lattice, it scales linearly with R. Thus 
we find that the heat of sublimation, which is already 
unusually high for buckyballs, increases even more with 
radius R. Using r — 4.56 and 77 = 2.05 for buckyballs 
yields AH — 58.46 kcal/mol; given that we use the near- 
est neighbor and the Derjaguin approximation (which is 
valid for 77 ^ 1), the agreement with the experimental 
value of 40.1 kcal/mol |2J] is surprisingly good. 
The bulk modulus follows as |E6i 
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For 77 ^ 1, Sq w 1 and B becomes independent of R. 
Therefore its scale must be set by e/a^ from dimen- 
sional considerations. The volume compressibility fol- 
lows as /3 = \/B. Extrapolating this to buckyballs, 
we find /3 = 3.65 x 10""'^^ cm^/dyn, which again is in 
surprisingly good agreement with the experimental value 
/3 = 7x 10-12 cmVdyn ||. 

The energy per area needed to cleave the crystal into 
two along lattice planes perpendicular to a given direc- 
tion is twice the surface energy 7 for this direction. For 
simplicity we will consider only the (lll)-direction. Be- 
fore doing so, we first discuss 7 for the layered material. 
Using the nearest neighbor approximation, we have to 
determine the energy per area at the equilibrium sepa- 
ration of two planar sheets. Taking the limit /i — > in 



Eq. (M), using x — '''/^ ^^^ minimizing for distance z 
yields 
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Again the combination e/tr^ follows from dimensional 
considerations. For carbon we find 7 = 150 erg/cm^, 
which compares quite well with the experimental value 
of 7 = 130 erg/cm^ for graphite [E3[. 

We now consider the (lll)-direction of the fcc-solid. 
Since a given molecule in a corresponding lattice plane 
has three nearest neighbors in the adjacent plane, 7 can 
be estimated to be 3Vb/2 times the area density in a 
close-packed plane with nearest neighbor distance r„„, 
2/V3r 
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For ?7 3> 1, So ~ 1 Ei-nd the surface energy 7 scales in- 
versely with R. Extrapolating to buckyballs, we find 
7 = 115 erg/cm^, that is a value smaller than the surface 
energy 70 of graphite. Although the interaction becomes 
stronger as i?, the number of vdW contacts decreases as 
1/R^ and the geometrical effects dominates. With in- 
creasing i?, it hence becomes easier to cleave the crystal. 
Note however that this does not lead to an enhanced 
roughening of the crystal surface. In fact the roughen- 
ing temperature scales as 70^ ~ R with a ^ R being 
the lattice constant |22]. Since the melting temperature 
scales with the potential depth which in turn scales with 
R, both temperatures scale with R and the roughening 
temperature remains a certain fraction of the melting 
temperature with varying R, as it usually is. 



IV. PHASE BEHAVIOR 

For the simple Lennard- Jones interaction of Eq. (Q), 
we can rescale the potential with the well depth, e, and 
the distance r with the hard core diameter a^ to write 
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which has no other free parameters. The resulting phase 
diagram in reduced temperature t — kT/e and vol- 
ume density $ = irpa^ /& (where p is density and a 
hard core diameter) is universal; the phase diagrams for 
all Lennard-Jones systems are identical when plotted in 
units oft and $ [law of corresponding states). In Sec. [I B 
we derived effective interaction potentials for single- and 
multi-walled nanoparticles. Single-walled nanoparticles 
interact via the Girifalco potential Vq from Eq. (g). Al- 
though derived from the Lennard-Jones potential from 
Eq. dl^) for two single atoms, it involves an additional 



length scale, the particle diameter 2i? = rja. Multi- 
walled nanoparticles interact with the complicated po- 
tential Vf given by Eqs. (|j,^). Here another length scale 
enters, the wall thickness h — va. Since they involve 
more than one length scale, the potentials Vq and Vp do 
not lead to a law of corresponding states and their phase 
behavior need s closer inspection. 

In Sec. [IC we showed that for large _R, the Derjaguin 
approximations given in Eq. (HO) and Eq. dl3) for Vg and 
Vp, respectively, give good approximations of the poten- 
tial well. In Sec. Ill this model was used to predict certain 



physical properties of solids which are essentially deter- 
mined by the properties of the potential near its mini- 
mum. In regard to phase behavior, the situation is more 
complicated, since now the entire range of the potential 
matters. However, recent work P,n,p[ suggests that for 
an isotropic interaction potential with hard core repul- 
sion and some attractive component for which potential 
depth and hard core diameter have been used to scale 
temperature and density to reduced units, the topology of 
the phase diagram will be essentially determined by one 
additional feature of the effective potential, namely its 
interaction range. This quantity is well defined only for 
square well potentials where it is the ratio of well width 
to hard core diameter. For continuouspotentials, the fol- 
lowing definition has been employed [[7|: 5 = {ri —ro)/rQ 
where tq is the position of the potential minimum and ri 
the distance at which the potential has fallen to 1/100 
of the potential depth at tq. Note that the interaction 
range S is independent of the potential depth e. For the 
Lennard-Jones potential given in Eq. (p9|), S — 1.42, and 
for the Girifalco potential for buckyballs given in Eq. (||), 
S = 0.83. In the framework of a variational scheme for 
the Double- Yukawa potential, it was found that the gas- 
liquid coexistence disappears for S < 0.4, and that an 
isostructural solid-solid coexistence appears for 6 < 0.02 
||7|. Similar values where found for the square well po- 
tential both in a simple van der Waals theory Q and in 
extensive MC simulations 0. 

The physical reason for the disappearance of the gas- 
liquid coexistence is well known. The fluid-solid phase 
transition is driven mainly by changes in free volume: 
at high densities, the solid becomes more favorable since 
its ordered structure provides more space for fluctuations 
(therefore hard sphere systems exhibit a entropy-driven 
fluid-solid transition). In contrast, the gas- liquid coexis- 
tence is driven by the attractive interaction: at low tem- 
peratures, the fluid phase separates so that the gas and 
liquid phases can profit entropically and energetically, 
respectively. The critical temperature for the fluid-fluid 
coexistence scales linearly with the depth of the poten- 
tial and also depends on its higher moments. Even if the 
potential is only slightly attractive, the critical point will 
exist by itself. However, the corresponding gas-liquid co- 
existence will survive in the overall phase diagram only 
if the critical temperature lies above the temperature at 
which the fluid coexists with the solid at the critical den- 
sity. In Fig. II we schematically depict the two possible 



outcomes which we have to consider in this work. 

The calculation of phase diagrams from various inter- 
action potentials has reached a high level of sophistica- 
tion and different schemes are readily available for a de- 
tailed analysis (compare [n3[ for an overview of different 
methods applied in the case of buckyballs). In the fol- 
lowing we ask how does the gas-liquid coexistence disap- 
pear from the phase diagram of hollow nanoparticles as 
a function of R and h. For this purpose, it is sufficient to 
adopt the simple criterion which was derived within the 
framework of a sophisticated variational scheme for the 
Double- Yukawa potential: the gas-liquid coexistence dis- 
appears if the interaction range S falls below S « 0.4. The 
basic justification for our approach lies in the fact that 
the interaction potentials under consideration can be fit 
well by the Double- Yukawa potential. This potential has 
in fact three parameters , one of which scales the po- 
tential depth. The remaining two can be used to adjust 
the position of the potential minimum and the interaction 
range. The zero crossing of the potential can be identified 
with the hard core diameter and for the Double- Yukawa 
potential equals one by definition. For colloids as well as 
for the hollow nanoparticles discussed in this work, the 
position of the potential minimum should always be close 
to the hard core diameter, that is little larger than one. 
Therefore we are left with only one degree of freedom, 
which in our approach can be used to fit the interaction 
range 6. The interaction range for the full potentials Vq 
and Vf can be determined numerically. The correspond- 
ing Derjaguin approximations can be used to gain some 
physical insight into the numerical results. In fact, for 
very large R the Derjaguin approximation should cor- 
rectly describe not only the part of the potential which 
contains the potential minimum, but also the part which 
contains the interaction range. 

We first discuss single-walled nanoparticles. The Der- 
jaguin approximation to the Girifalco potential Vg is 
given in Eq. dlQ). Combining it with the result for its 
potential depth, Eq. ([ll|), gives 
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Thus we find that for large R, the particle diameter, 
2R = r]a, simply serves to rescale the separation distance. 
With So from Eq. n^ and si — 1 + 4.56/?7 determined 
numerically, we find 
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Thus we conclude that the dimensionless interaction 
range S scales inversely with R for large R. This re- 
sult can be understood as follows: on an absolute scale, 
the interaction range is essentially constant since it corre- 
sponds to the range of the Lennard-Jones interaction be- 
tween atoms. It is made dimensionless by rescaling with 
the position of the effective potential minimum, which is 



fixed by the Lennard-Jones potential close to 2i?; there- 
fore the dimensionless interaction range S scales inversely 
with R. In Fig. O we plot the interaction range of the 
Girifalco potential and the result from the Derjaguin ap- 
proximation, Eq. (p^). Although the curves collapse onto 
each other for large 77, for rj — 2.05 they disagree consid- 
erably {S = 0.83 versus 6 = 1.27). Guerin (2^ fit the Gir- 
ifalco potential to a Double- Yukawa potential for which 
he found 6 = 0.56. The difference between this value for 
the fit and S = 0.83 for the Girifalco potential arises since 
the fitting procedure does not preserve the value of the in- 
teraction range. The approach used in this work amounts 
to using the three free parameters of the Double- Yukawa 
potential to fit it to the potential depth e, the position of 
the minimum vq and the interaction range S of the po- 
tential under investigation, thus preserving the value of 
the interaction range. Fig. |g shows that for 77 w 7, the 
full potential has S « 0.4 and the liquid-gas coexistence 
is expected to disappear. This corresponds to i? « 12 A, 
i.e. Cggo- On the basis of these results, we have to con- 
clude that buckyballs Ceo are in fact rather far away from 
loosing their gas-liquid coexistence. 

We now proceed to discuss multi- walled nanoparticles. 
The Derjaguin approximation is given in Eq. (0). Since 
it involves one more length scale, it cannot be treated 
in the same way as the Derjaguin approximation for the 
Girifalco potential. However, we showed above that the 
potential saturates for large h. The limit h ^ a is given 
in Eq. (p2[). Combining this with the result for the po- 
tential depth, Eq. Ifi3i), gives 
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As in the case for single- walled nanoparticles, 2R ~ yo' 
rescales the particles separation. With so from Eq. (|l^ ) 
and si = 1 + 66.1 8/77 determined numerically, we now 
find 

Thus the dimensionless interaction range 5 scales in- 
versely with R in both cases, although it has a larger 
prefactor in the multi- walled case, so that the gas-liquid 
coexistence disappears at larger values of R for larger 
wall thickness h. 

In Fig. M we plot the numerical results for the inter- 
action range S of the full potentials and their Derjaguin 
approximations. Different values of the interaction range 
S are shown as isolines in the [v, 77)-plane, and the critical 
interaction range Sc — 0.4 is drawn as a dashed line. The 
gas-liquid coexistence exists only in the region below the 
dashed line. Both diagrams agree with the general trends 
which were explained above: the coexistence disappears 
for large R and small h, and the critical line levels off 
for large h. However, the Derjaguin approximation over- 
estimates the interaction range, the real region with co- 
existence being in fact much smaller. From the results 



for the full potentials, we can conclude that the coexis- 
tence disappears for 77 « 7 (i? « 12 A) for single-walled 
and for r/ w 20 (i? w 35 A) for filled nanoparticles, re- 
spectively; the crossover with increasing wall thickness is 
rapid and on an atomic length scale. In particular we now 
can predict that buckyballs (single- walled, R = 3.55 A, 
•q — 2.05, h — V — Q) will have a gas-liquid coexistence, 
while typical metal dichalcogenides nanoparticles (multi- 
walled, i? = 600 1, ?7 = 346, /i = 93 1, j^ = 27) will 
almost certainly have none. 



V. ELASTIC DEFORMATIONS 

We now briefly discuss at which values of R and h 
the effects investigated above may be modified by the 
deformability of the particles. Although fullerene-like 
material interacts only through vdW forces, the result- 
ing elastic deformations can be substantial, as has been 
shown both experimentally and theoretically for carbon 
nanotubcs |28] . A detailed treatment of the elastic defor- 
mation of hollow nanoparticles is difhcult and relatively 
unexplored. Here we adopt a simple scaling approach 
which has been used before to investigate the mechanical 
stability of hollow nanoparticles in tribological applica- 
tions 1^. 

We ask when will considerable elastic deformations of 
hollow nanoparticles occur due to vdW interactions as a 
function of R and h. We consider two hollow nanopar- 
ticles adhering to each other due to their mutual vdW 
interaction as depicted in Fig. ^. The hollow nanoparti- 
cle is assumed to be an elastic shell with prefered radius 
R. Its deformation energy has two parts: bending energy 
characterized by the bending constant n and stretching 
energy characterized by the two-dimensional Young mod- 
ulus G. If an external force is applied, the balance be- 
tween these two contributions leads to a localization of 
the deformation [^. For small adhesive load, the shell 
flattens in a contact region and the overall deformation 
energy scales as 



Edef ^ H 



(24) 



where H is the indentation (compare Fig. ph. The vdW 
energy gained on deformation can be estimated as fol- 
lows: apart from some numerical prefactor, the adhesion 
energy per area is essentially ^tja^ (this follows from 
minimizing for distance z in Eq. (|7|); the limit of single 
walls has been given in Eq. (O)). For small _ff, it fol- 
lows from simple geometry that the area over which the 
two shells approach each other scales with RH (compare 
Fig. H) . Therefore we have 



£^, 



vdW 



x'Re 



H 



(25) 



For single-walled nanoparticles, x ti3,s to be replaced by 
T, which however has nearly the same value. We now can 



solve Edef = Eydw for indentation H. We then ask for 
the critical radius Re at which the indentation becomes 
considerable, i.e. of the order of the smallest length scale 
in our problem, which is a. This yields 
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For single-walled fuUerenes, the value for the bending 
rigidity k can be extracted from molecular calculations to 
be K = 1.6 X 10~^^ erg. The value for the two-dimensional 
Young modulus G can be extracted from the elastic mod- 
uli of graphite as G = 3.6 x 10^ erg/cm^. Then we find 
rjc = 16.5, that is considerable deformation is expected 
for R > 30 A. For single- walled nanoparticles made from 
metal dichalcogenides M0S2 and WS2, G is smaller by a 
factor 4 and k is larger by a factor 10, thus ijc increases 
slightly (by a factor 1.3). 

For nanoparticles with a few walls {a < h < R), one 
can show in the framework of linear elasticity theory that 
K = Ciih^/12 and G = Cnh where Cn is the in-plane 
stretchin g el astic constant of the corresponding layered 
material |^^. Using these relations yields 
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Thus the critical radius Re scales linearly with the wall 
thickness h. The values for Cn are 1060, 238 and 150 x 
W° erg/cm3 for C, M0S2 and WS2, respectively [||. 
For WS2 particles with h = 100 A, the critical radius 
following from Eq. ( pTJ ) is well above 1000 A. 

We found above that the gas-liquid coexistence disap- 
pears for Rk, 12 a for single- walled and for i? « 35 A for 
filled nanoparticles. Our estimates for the onset of elastic 
effects is at least of the same order of magnitude. There- 
fore we conclude that the elastic effects will not affect hol- 
low nanoparticles which do have a fluid-fluid coexistence. 
In general, the effect of wall thickness h of multi-walled 
nanoparticles is much stronger on the elastic response 
than it is on the phase behavior. Therefore multi-walled 
variants are favorable if one wants to make use of the 
unusual properites of hollow nanoparticles without being 
restricted by elastic effects. 



VI. CONCLUSION 

In this work, we predicted the phase behavior and ma- 
terial properties of hollow nanoparticles as a function of 
radius R and wall thickness h. The synthesis of hol- 
low nanoparticles is a rapidly developing field, driven by 
the special properties of particles in the nanometer range 
which may allow many new applications to be developed 
in the future flj] . Although our analysis is aimed at hol- 
low nanoparticles made from anisotropic layered material 
like carbon (fuUerenes) or metal dichalcogenides (inor- 
ganic fuUerenes), it could also be applied to colloidal or 



biological examples of hollow nanoparticles. Our starting 
point was the Lennard-Jones interaction between single 
atoms, which we integrated over the appropriate geome- 
tries in order to obtain effective potentials for the vdW 
interaction between single- and multi-walled nanoparti- 
cles. The subsequent analysis was twofold: Derjaguin 
approximations for the effective potentials provided the 
correct scaling laws and physical insight, while numeri- 
cal investigations of the full potentials allowed us to make 
quantitative predictions. 

We first showed that crystals of hollow nanoparticles 
have unusual material properties. Since the effective con- 
tact area and therefore the potential depth scales linearly 
with radius R for large R, their heat of sublimation and 
surface energy will scale linearly and inversely with R, 
respectively. This means that with increasing R, it will 
become more difficult to melt, but easier to cleave the 
crystal. The thickness h has little influence here since the 
vdW interaction saturates over the atomic length scale 
of the Lennard-Jones potential under close-packed con- 
ditions. For reasonable values of R, we predicted values 
for the heat of sublimation and surface energy which for 
vdW solids are unusually high and low, respectively; this 
was already found experimentally for fuUerite, but should 
be more pronounced for larger fullerenes and inorganic 
fuUerenes. 

Our discussion of phase behavior centered around the 
existence of a liquid phase and the concept of the inter- 
action range. We argued that after scaling out hard core 
diameter and potential depth, the interaction range is 
the relevant feature of an effective potential which deter- 
mines phase behavior. We showed that for large radius R, 
the interaction range scales inversely with R, and again 
is little affected by thickness h. The numerical analysis 
supported this finding and allowed to calculate the inter- 
action range as a function of R and h. Using recent work 
on the Double- Yukawa potential, which can be fit well 
to our potentials, we could identify the disappearance of 
the liquid phase with a dimensionless interaction range 
S « 0.4. We then found that the liquid phase will disap- 
pear for i? « 12 A for single- walled and for i? « 35 A for 
filled nanoparticles. Although several theoretical studies 
investigated the phase behavior of buckyballs as a likely 
candidate for the non-existence of the liquid phase in a 
non-colloidal system, we conclude that buckyballs with 
R = 3.55 A are in fact rather far away from losing their 
liquid phase. 

Finally, we showed that the wall thickness h has a 
strong effect on the elastic properties of hollow nanopar- 
ticles. Using scaling arguments for elastic shells with 
prefered radius R and the known scaling of the elastic 
moduli with h, we predicted that elastic effects should 
not interfere with the phase behavior of those (inorganic) 
fullerenes which do have a liquid phase. 

Despite the widespread theoretical interest in the 
phase behavior of buckyballs, little experimental work 
has been done on the phase behavior of hollow nanopar- 
ticles. Since the potential depth scales linearly with the 



radius R, the corresponding energy scales are large and 
experiments have to be conducted at temperatures of sev- 
eral thousand Kelvin, where their mechanical stability 
becomes a limiting factor. However, the temperature 
scale can in principle be lowered by up to two orders 
of magnitude by dispersing the nanoparticles in a fluid 
medium of high dielectric constant |2^ . Regarding hol- 
low nanoparticles which are impermeable to solvent, it 
would not be difficult to extend our analysis to the case 
of vdW interactions between spatial regions with three 
instead of two different dielectric constants. Our analy- 
sis should be directly applicable to the vdW interaction 
in systems with hollow polyelectrolyte and silica shells, 
whose walls are permeable to solvent [0,nS]. However, in 
this case additional effects like electrostatic interactions 
might have to be taken into account. Moreover, elastic 
and even plastic effects (low yield strength due to non- 
covalent bonding in the shell) are expected to be more 
prominent for these systems. 
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APPENDIX A: INTEGRATION OF 

LENNARD-JONES POTENTIAL OVER TWO 

BALLS OF UNEQUAL RADII 

In order to integrate the Lennard-Jones potential over 
two balls of unequal radii i?i and i?2 and center-to-center 
distance r, it is helpful to adopt the following coordinate 
system: consider a sphere of radius R centered at z ^ 
r on the z-axis. This sphere can be decomposed into 
caps of spherical shells of radius r' centered around the 
origin. Integrating out the two angles on the spherical 
caps leaves us with the following integral (20| : 



dU/(r') = - 
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E.g. for /(r) = 1 we obtain the sphere volume AttR^ /4. 
To integrate a given /(r) over two spheres, the same inte- 
gral has to be used twice, where in the second integration 
/(r) is the result of the first integration. For f{r) — l/r^ 
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FIG. 1. Schematic drawing of two multi-walled hollow 
nanoparticles with center-to-center distance r, particle radius 
R and wall thickness h. 



V 



0.4 

0.0 

-0.4 

-0.8 











■ 




/ 


h 
1 1 
1 1 

ft 

ii 
fi 
ji 

r ) 
1 ( 
t r 


1 / 
I ' 
' / 



0.5 



1.5 



0.05 



V[kT] 



-0.05 



-0.1 



V[kT] [ 





-10 








(b) 




1 / 

* / ' 



10 



12 



14 

r[A] 



FIG. 2. Different van der Waals interaction potentials: the 
Lennard- Jones potential for carbon atoms (a) is two orders 
of magnitude weaker than the Girifalco potential for buck- 
yballs (b, solid). A Lennard- Jones potential with the same 
zero-crossing and potential depth as the Girifalco poten- 
tial (b, dashed) shows that the Girifalco potential is rather 
short-ranged. The vertical line indicates that the potential 
diverges at a finite value of separation. 



FIG. 3. Interaction potential for single-walled nanoparti- 
cles: Girifalco potential (full) and Derjaguin approximation 
(dashed) for different values of 77 = 2R/a: from right to left 
ri = 2.05 (buckyballs) , 4 and 10. All potentials are scaled 
such that the Derjaguin approximation has its minimum at 
— 1. With increasing rj, the potential well of the full potential 
falls into the range described by the Derjaguin approximation. 
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FIG. 4. Interaction potential for multi- walled nanopar- 
ticles: Full potential (solid) and Derjaguin approximation 
(dashed) for rj — 2R/a = 10 and f — h/a — 0.8 (lower 
curves) and 0.4 (upper curves). All potentials are scaled with 
the same factor as in Fig. H. The Derjaguin approximation 
improves with increasing wall thickness v. 
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FIG. 6. Double logarithmic plot of the interaction range 5 
of single-walled nanoparticles as a function of r; = 2R/a for 
the Girifalco potential (full) and its Derjaguin approximation 
(dashed). The dotted lines at 5 = 0.4 and S = 0.02 separate 
the long-, intermediate- and short-ranged regimes. Bucky- 
balls with rj — 2.05 still have the gas-liquid coexistence which 
disappears only for S < 0.4, that is 77 > 7 (i? > 12 A). 



FIG. 5. Phase diagrams (a) with and (b) without gas-liquid 
coexistence as a function of reduced temperature t = kT/e 
and volume fraction 4> ~ 7i"po"'^/6. G, L, S and F denote 
gaseous, liquid, solid and fluid phases, respectively. In (b) 
the attractive interaction potential is so short-ranged that 
the gas-liquid coexistence (although existing by itself) is sup- 
pressed in the overall phase diagram by the fluid-solid coexis- 
tence. For a Double- Yukawa potential, the crossover between 
the two topologies of the phase diagram occurs when the in- 
teraction range decreases below the critical value 5 ~ 0.4. 



12 



100 



11 



80 



60 



40 



20 





FIG. 8. Schematic drawing of two multi-walled hollow 
nanoparticles deformed due to vdW adhesion. H is the inden- 
tation. It follows geometrically that the contact area scales 
with RH. 
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FIG. 7. Interaction range 5 of multi-walled nanoparticles 
as a function of 77 = 2R/a and v = h/a for (a) Derjaguin ap- 
proximation and (b) full potential. The dotted line delimits 
the allowed region h < R, that is r; > 2^. The dashed line 
marks d = 0.4. Only the region between the dotted and the 
dashed lines has a gas-liquid coexistence, whose size is overes- 
timated in the Derjaguin approximation. Filled nanoparticles 
(77 = 2u) do not have a liquid phase for rj > 20 {R> 35 A). 
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